A second-order phase transition in CDT 
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Abstract 

Causal Dynamical Triangulations (CDT) are a concrete attempt to define a non- 
perturbative path integral for quantum gravity. We present strong evidence that 
the lattice theory has a second-order phase transition line, which can potentially 
be used to define a continuum limit in the conventional sense of nongravitational 
lattice theories. 
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1 Introduction 



Dynamical triangulations (DT) were invented as a nonperturbative regularization 
of bosonic string theory and thus also of two-dimensional quantum gravity coupled 
to conformal matter. This program was both a failure - in showing that even in a 
nonperturbative setting no bosonic string theory exists in dimension two or larger 
- and an amazing success, in providing a versatile regularization of 2d quantum 
gravity coupled to conformal matter with central charge c < 1, i.e. noncritical 
string theory. Surprisingly, in many ways the regularized theory turned out to be 
easier to solve analytically than the corresponding continuum theory. 

Encouraged by this, DT was generalized to provide a regularization of quan- 
tum gravity in three [UEIE] and four dimensions [UH]. The naive expectation was 
that if a 'stand-alone' four-dimensional theory of quantum gravity existed, the 
regularized theory should have a second-order phase transition, which could be 
used to define a continuum theory of quantum gravity. Second-order transitions 
are usually characterized by a divergent correlation length associated with prop- 
agating field degree (s) of freedom, in the case at hand presumably of a gravitonic 
nature. If a given phase transition point was a UV fixed point, one could also 
attempt to make contact with Weinberg's asymptotic-safety scenario, for which 
plenty of corroborating evidence has been found recently (see [6] for reviews). 

A phase transition point was indeed located in 4d DT, and at first believed to 
be of second order [HEHS]- However, analyzing larger lattice systems changed the 
verdict to a first-order phase transition where no obvious continuum limit could be 
defined, at least not when using the Regge version of the Einstein-Hilbert action 
[9l [10] . Neither did one find convincing evidence of a good classical behaviour of 
large-scale geometry away from the phase transition. 

Partly triggered by this impasse, a modified lattice model in terms of Causal 
Dynamical Triangulations (CDT) was proposed, and subsequently shown to have 
long-distance properties in agreement with (semi-) classical gravity [13J. It still 
uses the Regge-Einstein-Hilbert action, but assumes the existence of a global 
time-foliation, and has a more complicated phase diagram than the simplest DT 
mo del0 In what follows, we will provide strong, new evidence that - unlike its 
Euclidean counterpart - 4d CDT quantum gravity possesses a second-order phase 
transition line. 

2 Causal Dynamical Triangulations 

We begin with a brief account of CDT, focusing on several important aspects 
(see [TU [15j [16j [17] for technical details and [18] for reviews). CDT can be char- 
acterized as a nonperturbative path integral which is as close as possible to a 
canonical quantization: spacetime histories share a global foliation, where each 

1 The notable similarities with 'Hofava-Lifshitz gravity' [TT| are being explored, see e.g. [12]. 



2 



leaf is a spatial hypersurface, given in terms of a three-dimensional triangulation 
of fixed topology T, built from equilateral tetrahedra with link length a s , and 
labeled by a discrete proper time t n . Adjacent hypersurfaces are connected by 
four-simplices, resulting in spacetime histories of the form of four-dimensional tri- 
angulations of topology T x [0, 1]. We use T = S 3 and impose periodic boundary 
conditions in time, such that the spacetime topology is S 3 x S 1 . The geometry 
of each spacetime is fixed by how the four-simplices are 'glued together' to form 
a simplicial manifold, and by the lengths of its links, which come in two types: 
spacelike links which lie entirely within a given hypersurface, and timelike links 
whose endpoints lie on adjacent hypersurfaces, and which have a (squared) edge 
length a\ = aa 2 s) for some fixed relative scaling parameter a < 0. 

The CDT gravitational path integral is a sum over all geometrically inequiva- 
lent triangulations of this type with a fixed number of time steps, with amplitudes 
depending on the above-mentioned Regge-Einstein-Hilbert action [19] . Since in 
four dimensions analytical methods are mostly unavailable we use Monte Carlo 
simulations to extract physical results. To do this we must convert the path in- 
tegral to a statistical partition function by applying a Wick rotation, which due 
to the foliated structure exists globally [T5]. It can be implemented by rotating 
a — > —a in the lower-half complex plane, leading to the Euclidean Regge action 

SE = hI + 2A) ^ + 6A)iVo + KiNi + A(iV4 + ^ 4,1)) ' (1) 

where Nq and N± denote the total number of vertices and four-simplices and ivj 4 ' 1 ^ 
counts the subset of four-simplices which have four vertices on one hypersurface 
and the fifth one on a neighbouring one. The couplings k,q, k a and A are linearly 
related to the bare inverse gravitational coupling, the bare cosmological coupling 
and the parameter a introduced above. Redefining R A = k a + A, we obtain the 
Euclidean Regge action implemented in the computer simulations, namely, 

S R egge = -(«o + 6A)iV + k 4 N A + AAr| 4,1) =: -k N + k A N A + A conj(A), (2) 

where for later convenience we have introduced the quantity conj(A) = N^' 1 ^ — 
6Nq conjugate to A. This turns the gravitational path integral into a statistical 
partition function with Boltzmann weights exp(— S Regge)- We will use the freedom 
to switch to a different ensemble, obtained by keeping N A (measuring the four- 
volume of the system) fixed, instead of its conjugate k a . In this way we can treat 
A^4 as a finite-size scaling parameter which does not appear in the phase diagram 
of the putative continuum theory. The remaining couplings k,q and A span the 
phase diagram, which we will go on to explore in the next section. 

3 The phase diagram of CDT 

A qualitative description of the CDT phase diagram first appeared in [TJ], with 
a quantitative plot presented later in [12]. The diagram exhibits three phases, 
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Figure 1: The phase diagram of CDT. The large crosses represent actual mea- 
surements. 

labeled A, B and C in Fig. [TJ A geometric characterization of the phases can 
be given in terms of their distinct spatial volume profiles N 3 (t), measuring the 
three- volume in lattice units as a function of proper time t. As described in detail 
elsewhere [TO [T6] , the average large-scale geometry found in phase C shows the 
scaling behaviour of a genuine four- dimensional universe. The average volume 
profile matches beautifully that of a Euclidean de Sitter spacetime. Even the 
quantum fluctuations around this emergent background geometry are described 
well by a cosmological minisuperspace action [T7J [20] • The situation in the other 
phases is completely different. The typical volume profile of a configuration in 
phase A shows an almost uncorrelated sequence of spatial slices, while the con- 
figurations in phase B are characterized by an almost vanishing time extension. 

A preliminary analysis in [12] suggested that the A-C transition is of first 
order, similar in nature to the first-order transition observed in the DT formalism 
mentioned above. However, as also pointed out in [21], to nail down this result 
the numerical evidence still needs to be strengthened. By contrast, our interest in 
the present letter is in the order of the B-C transition, about which considerable 
doubt remained in [12] . We will present strong evidence below that it is a second- 
order transition line. 

Before doing so, it is instructive to analyze how the transitions change as 
we move along the respective phase transition lines, while holding the system 
size fixed. The A-C transition line is characterized by a jump in Nq, the variable 
conjugate to the coupling constant k,q, when we cross the line by changing Kq- We 
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see no appreciable change in this signal when we move along the A-C transition 
line, although we have not examined closely the triple point where all three phases 
meet. 

The situation is very different for the B-C transition. When changing the 
coupling constant A and crossing the transition line, we observe a jump in the 
variable conj(A). However, moving to smaller values of k on the left, the jump 
decreases. Around k,q = 1.0, no signature of a phase transition remains. We 
conclude that the B-C transition has an endpoint, which for N4 = 80k is located 
around k,q = 1.0. Moving to the right the jump increases, and around k,q = 2.3 
the transition is so strong that we get stuck in metastable states. The dashed 
part of the B-C line in Fig. [T] marks the region where conventional methods 
are insufficient to measure the location of the phase transition with acceptable 
accuracy. We are analyzing currently whether the use of multicanonical Monte 
Carlo simulations can help in resolving this issue. 



4 The order of the B-C phase transition 

Measuring the order of phase transitions requires some care. To confirm that a 
phase transition is not a first-, but a second-order transition, one can try to mea- 
sure various so-called critical exponents. One such exponent measures the shift of 
a transition point with system size. Recall first how this works for a conventional 
system such as the Ising model with volume V = L d , where d is the dimension 
of the system [22]. Considering the Ising model's temperature-driven phase tran- 
sition and using the location of the maximum of the magnetic susceptibility to 
define a transition point j3 c (V), one finds the power-law behaviour 

|/3 c (oo) -P C {V)\ oc V- 1/ud (3) 

for sufficiently large system sizes. The exponent v governs the increase of the 
correlation length in a second-order transition as one moves towards the critical 
point /3 c (oo) on an infinite lattice. For first-order transitions there is no correlation 
length and one expects a scaling like ([3]), with vd replaced by an exponent v where 
v — 1 [23]. A sufficiently strong violation of v = 1 therefore signals the presence 
of a second-order transition. Another quantity of interest is the Binder cumulant 

associated with an observable O [23] , which is always nonpositive and zero if the 
histogram of O is Gaussian. Evaluating B as a function of the couplings, its 
local minima lie at transition points. We can measure these minima for different 
system sizes and by extrapolation determine B^il/N^ = 0). At a second-order 
transition the histogram of O should converge to a single Gaussian distribution 
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Figure 2: Measuring the location A c of B-C transition points at k = 2.2 for 
different system sizes A/4 to determine the shift exponent v. 

with the Binder cumulant going to zero. At a first-order transition it will go to a 
nonzero constant, however, its value at a weak first-order transition can be small. 

To analyze the B-C transition, we fixed k,q = 2.2 and used systems of size 40, 
50, 60, 80, 100, 120, 140 and 160k. The number of sweeps used was approximately 
2.5 • 10 6 , with one sweep consisting of one million attempted Monte Carlo moves. 
We have measured the shift exponent v for the asymmetry parameter A using 

A c (iV 4 ) = A c (oo) - CN^'\ (5) 

where C is a proportionality factor, and A c has been defined using the location 
of the maximum of the susceptibility Xconj(A) = (conj(A) 2 ) — (conj(A)) 2 . Fig. [2] 
shows the measured data points (error bars too small to be included) and the 
best fit through all of them, yielding v = 2.39(3). To judge whether our range 
of system sizes lies inside the scaling region we have performed a sequence of 
fits by successively removing the data points of lowest four-volume, leading to 
//-values 2.51(3), 2.49(3) and 2.51(5), where the last fit was performed with all 
but five data points removed. This suggests that the data point with the lowest 
four-volume lies outside the scaling region. Removing it from the fit we obtain 

£ = 2.51(3). (6) 
This result makes a strong case for a second-order transition, since the prediction 
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Figure 3: Dependence of the minimum of the Binder cumulant -B C onj(A) on the 
(inverse) system size at the B-C transition. At a second-order transition, B mm — > 
in the infinite- volume limit. (Fit excludes the two points on the right.) 

v — 1 for a first-order transition is clearly violated. (By contrast, for the A-C 
transition one finds v w 1 as will be reported elsewhere [21].) 

Lastly, we have investigated how the minimum of the Binder cumulant 
depends on the system size. Fig. [3] shows -B^n^A) as a function of inverse system 
size (errors approximately equal to the dot radii). Inside the scaling region the 
minimum of the cumulant is expected to have a power-law behaviour. To under- 
stand which data points lie inside the scaling region, we have again performed a 
sequence of fits by successively removing the points of lowest four-volume. This 
has led to the exclusion of the data points for N± = 40k and iV 4 = 50k from the 
fit shown in Fig. [3J 

Table [T] collects the extrapolations Bq 1u (N4 — > 00) for several observables O. 
As indicated earlier, it is always difficult to make a strong case for a second- 
order transition based on Binder cumulant measurements alone, because weak 
first-order transitions may show a convergence to a nonzero value close to zero. 
Nevertheless, in the case at hand all our measurements are mutually consistent 
and in excellent agreement with the limiting value 0, further corroborating our 
claim of the second-order nature of the transition. 
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observable O 


B™ m (N A ->■ oo) 


conj(A) 


-0.003(4) 




-0.001(3) 


N 2 


-0.0000001(3) 




-0.000003(7) 


N 


0.0000(3) 



Table 1: Measurements of B^^N^ — > oo) for various observables O, where Nk 
denotes the number of /c-dimensional (sub-)simplices in the triangulation. 

5 Discussion 

We have succeeded in our goal of determining the order of the B-C transition in 
CDT quantum gravity by applying two distinct methods, namely, measuring the 
shift exponent and analyzing Binder cumulants. The measured shift exponent 
u = 2.51(3) represents a strong violation of the prediction v = 1 for a first- 
order transition. Also the results of the Binder cumulant analysis are clearly and 
unambiguously consistent with the second-order nature of the transition. 

From this we conclude that there is strong evidence that the B-C transition is 
of second order, making four-dimensional CDT quantum gravity the first known 
instance of a dynamically triangulated model (without matter coupling) in any 
signature and dimension which displays such a transition. This result is poten- 
tially very attractive. It opens the door to studying critical phenomena in CDT 
and defining a continuum limit where the lattice spacing (the UV cut-off) is 
taken to zero, just as one does in standard lattice quantum field theories with 
nondynamical geometry. 
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